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Abstract 

We calculate in this article the transport coefficients which characterize the 
dynamics of solitons in quantum field theory using the methods of dissipative 
quantum systems. We show how the damping and diffusion coefficients of 
soliton-like excitations can be calculated using the integral functional formal- 
ism. The model obtained in this article has new features which cannot be 
obtained in the standard models of dissipation in quantum mechanics. 

PACS numbers: ll.lO.-z, ll.lO.Lm, 03.70.+k, 05.40.+j, 05.60.+W, 03.65.Db 
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I. INTRODUCTION 



There are still controversies about the calculation of transport properties involving 
soliton-like particles (see and references therein). It is not completely clear how the 
thermalization of collective excitations, such as solitons, can be carried out. In general, 
the pathway is based on heuristic approaches which take into account averages over con- 
figurations without going deeply inside the dynamical characteristics of the formation of a 
collective excitation. 

The main goal of this article is to link two fields of study which seem to be far apart up 
to now but, as we shall show, are really closely related. 

One of these areas is the quantum theory of solitons developed in the seventies which 
teaches us how to quantize classical theories which have solitons as solutions of their equa- 
tions of motion. The other area is the field of dissipative quantum systems which had a 
great development during the eighties |^. 

When we study the dynamics of a classical particle moving through a viscous environment 
we learn that one can describe its motion using two parameters, the damping and the 
diffusion coefficient. The former is related to the systematic force applied on the particle by 
the environment and the latter is related with the fluctuations due to the interaction (for a 
clear discussion about this issue see |^). 

It was shown some years ago that the same physical situation occurs in quantum systems 
P], therefore if a particle interacts with an environment such as a set of decoupled harmonic 
oscillators (the oscillator model) its motion is damped, which means in quantum mechanical 
terms that the center of motion of its wavepacket undergoes a damped motion while its 
width increases with time. Exactly as in the classical problem all the transport properties 
of the quantum system can be described by a damping and a diffusion coefficient. 

One of the main results of the standard approach is that the damping time is independent 
of temperature and even at zero temperature the motion of the particle is damped. As 
will show in this article this behaviour is not expect for solitons since the soliton and the 
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environment have the same microscopic origin. 

From the point of view of quantum field theory, classical and translational invariant 
theories with solitons can be quantized by the collective coordinate method 0]. In the 
classical theory, solitons are field configurations which move in space without changing its 
shape. Therefore, the knowledge of the position of the center of soliton as a function of time 
is enough to describe its dynamics. In quantum field theory the center of the soliton plays a 
special role. The center of the soliton can be viewed as a true quantum dynamical variable 
and we describe the soliton as a particle. The problem here is that at finite temperatures not 
all the degrees of freedom (infinite, indeed) collaborated in the formation of the soliton and 
therefore this soliton is never free to move as in the classical theory, that is, there is always 
a residual interaction due to quantum degrees of freedom which changes drastically the 
dynamical properties of the system. Our task is to show that, under some approximations, 
this residual interaction gives rise to a damped motion of the soliton, in other words, the 
soliton undergoes a brownian motion. 

Form this point of view we are presenting here a new model for dissipation in quantum 
mechanics and, as we would expect, very different from the oscillator model described before. 

We could have started from the definition of the model , which by itself is completely 
new in the field of quantum dissipation, without mentioning its origin in quantum field 
theory, but we think it is of interest to show the generality of the model in this context. 

We will show here how we can calculate the transport properties in this framework, 
namely, the damping and the diffusion coefficient. We will show that the damping time 
depends on the temperature, a result which can not be obtained in the oscillator model. 

In section II we present the results of the quantization of the solitons in field theory and 
the approximations involved in our model. Although the same approach can be found in 
great detail in Ref. we include this topic here for the sake of completeness. In Section III 
we develop the the integral functional formalism which allows us to calculate the transport 
properties. Section IV contains our conclusions. 
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II. QUANTIZATION OF SOLITONS 



Let us consider a classical action in 1+1 dimensions for a scalar field, 

S^JcPx {d^<f>) {d^<t>) - U{g, 0)) (2.1) 

where C/ is a potential function of the field (f) and the coupling constant g. Here we will 
assume that U has a degenerate absolute minimum. We can always rescale the field such as 
4> ^ 4>/g and the potential as U U/g^ (by naive dimensional analysis) and rewrite (2.1) 
as 

S^^Jcfx^^ {d^cl^) (a^</)) - [/(</))) . (2.2) 
Consider the static part of (2.2), 

y = i!l'{''^\U^X + U(<l>)\ (2.3) 



2 \dx^ 

and its extremum value given by SV/S4>s = 0, namely, 

£(j)s ^ dU{(t>s) 
dx^ d(l)s 

which, due to the translational invariance of (2.1), must have the form 



(2.4) 



(ps = 4>s{x - xo) (2.5) 

where Xq is a constant. Eq.(2.5) is a static localized solution of the classical equations of 
motion. 

We expand (2.3) around the extremum, (2.4), in terms of the coupling constant g as, 

(l){x,t) ^ (l)s{x - xo) + g5(f){x,t) (2.6) 

and one finds, 

V[5cP] = y [0,] + I d'x S<P{x, t)(^-^+ j Scf>{x, t) (2.7) 

plus higher order terms in g. 
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Using (2.4) we get, 
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M^mi-^/<iV(^j (2.8) 

is the classical soliton mass. 

Therefore, up to leading order in g, the eigenmodes of the system are given by 

+ 1 1 *n(^ - ^O) = ^I'^nix - Xo). (2.9) 



Observe that the function d(j)s/dx is an eigenfunction with eigenvalue zero (use (2.4) in 
(2.9)) indicating that we should expect divergencies in high order terms in the perturbative 
expansion. The way to deal with this "zero" mode is the following: suppose we displace the 
center of the soliton, xq, by an infinitesimal amount Sx^. Thus, up to first order in Sxq, 

So<Ps = <Ps{xo + Sxq) - (psixo) = Tr-^5xo (2.10) 

OXq 

but since 0s is a function of the difference x — xo we can rewrite (2.10) as 

5o0, = -^5xo (2.11) 

therefore the eigenmode d(ps/dx is related with the movement of the localized solution, that 
is, the motion of the sohton. Call, 

where is a normalization constant. 
We rewrite 6(f) in (2.6) as, 

CXD 

S(j>{x, t) = J2 an{t)'^n{x - Xo). (2.13) 

n=0 



Therefore, from (2.6), 



4>{x,t) = 4>s{x - Xo) + j^ao{t)^ + gJ2an{t)'^n{x - Xo). (2.14) 

Notice that the first two terms in r.h.s. of (2.14) appear in the expansion obtained in 
(2.10). In this way we see that the center of the soliton is a true dynamical variable and we 
rewrite the expansion (2.14) as. 



t) = (f)s{x - xo{t)) +gY. cin{t)^{x - Xo{t)) (2.15) 

n=l 

which is known as the collective coordinate method. All the physics of this system is present 
in (2.15). The soliton, described by the motion of its center, is a collective excitation (since is 
a solution of the field equation) but it is coupled to all other modes by the relative coordinate 

X — Xq. 

We can now rewrite the action (2.2) in the presence of the soliton using (2.14) (see 
for details) and obtain the classical Hamiltonian for his problem, 

" ^ „£. ^ ' S (t + <^-^«) 

plus higher order terms in g. P is the momentum canonically conjugated to Xq, and Pn 
are the conjugated pairs for the modes and 

Gmn = ~ f dx^{x)^n{x) (2.17) 

J ax 

couples the soliton to the other modes (mesons) in the system. 

It is clear now that the soliton can not move freely, the mesons are scattered by the 
soliton which will move as a brownian particle. Although (2.16) is perturbative in the 
coupling constant, we will assume that the main physics of the problem is present in that 
expression. 

The quantization of (2.16) is straightforward. We impose the commutation relations 

[xo,P]=ih (2.18) 

[qn,Pm] = iflSnm (2.19) 

and write (2.16) in the symmetrized form (apart from the constant mass term) 

1 2 °° 

H = — (P-P^,) +Y.nnJX (2.20) 

^^'^ n=l 

where 
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which obeys the usual commutation relations for bosons, 

bn, hla = Sum ] K, &m = 0. (2.22) 

In (2.20), Pme is the momentum of the mesonic field which is given by 

Pme = Pd + Pnd (2.23) 



where 



P^= ^ (i?-) + (t^) Gmnblb^ (2.24) 



n,m=l 2^ \Vfim/ V 



Pnd= tA (tT-) - (tF) G^^n {hJn - hlhi) . (2.25) 



and 



The Hamiltonian (2.20) is well known in works on polaron dynamics and many 
perturbative approaches has been used in order to understand its physical contents. 
Observe that the Pnd does not commute with the total number of mesons 

oo 

N=Y. bik (2.26) 

n=l 

and therefore it breaks the number conservation. Actually this term is related with high 
frequency oscillations in which we are not interested. Our task is to show that in the sector 
of the Hilbert space where the number N is conserved the physics of (2.20) is that of a 
brownian particle. Therefore our new model in the context of quantum dissipation is based 
on the following Hamiltonian |^ 



2 

oo 



H=7kf-{P- E ^9mnblk ) : + E ^^nHk (2.27) 

n,m=l I n=l 



2M 

where : ... : means normal order and 



9mn = l:{ {^\" + {^\"\ G_ (2.28) 



is the coupling constant of the theory. 

We beheve that Hamiltonian (2.27) describes the main physics of the dynamics of sohtons 
at low energies and, by itself, can be used to model other physical systems where the coupling 
of the particle and the environment can not be described by the standard model. 

III. FUNCTIONAL INTEGRAL METHOD 

The starting point for the calculations of the transport properties of the soliton is the 
well-known Feynman- Vernon formalism that the authors have recently applied ^ to the 
Hamiltonian (2.27). 

We are interested only in the quantum statistical properties of the soliton and the mesons 
act only as a source of relaxation and diffusion processes. Consider the density operator for 
the system soliton plus mesons, p{t). This operator evolves in time according to 

p{t) = e-'^*/^p(0)e'^*/'' (3.1) 

where H is given by (2.27) and p{0) is the density operator at t = which we will assume 
to be decoupled as a product of the soliton density operator, ps{0), and the meson density 
operator, Pr{0), 

p{0) = psmdO) (3.2) 

where the symbol S_ refers to the soliton (system of interest) and R to the mesons (the 
reservoir of excitations). 

We will consider that the mesons are in thermal equilibrium at t = 0, that is, 

Pr{0) = (3.3) 

where 

Z = tvR (e-^^«) (3.4) 

with 
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Here trji denotes the trace over the mesons variables and Kb is the Bohzmann constant. 
Hb is the free meson Hamiltonian which is given by the last term on the right hand side of 
(2.27). 

As we said, we are interested only in the quantum dynamics of the system S, so, we 
define a reduced density operator 

Ps{t)=trn{p{t)) (3.6) 

which contains all the information about S_ when it is in contact with R. 

Projecting now (3.1) in the coordinate representation of the soliton system 

Xo\q>= q\q> (3.7) 

and in the coherent state representation for bosons (the mesons), 

bn I an >= On I a„ > (3.8) 

we have [0 (see ref. for details) 

Ps{x,y,t) = dx dy'J{x,y,t;x',y',0)p,{x',y',0) (3.9) 



Here we have used (3.2), (3.3), (3.6) and the completeness relation for the representations 
above, namely, 

Jdq\q><q\=l (3.10) 

/ \a><a\=\ 3.11 

J vr 

where (i^a = d{Rea)d{Ima) as usual. 

In (3.9), J is the superpropagator of the soliton, which can be written as 

J = r Dx r Dy e^(^°["]-^°t^l)F[x,y] (3.12) 

J x' Jy' 
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where 



■,A^iy\^\ (3.13, 



is the classical action for the free particle. F is the so-called influence functional, 

r A'^P, r ff^Ri _ _ |3|2 i3',2 



(3.14) 



where /? denotes the vector {Pi, P2, P3, Pn) and 5*/ is a complex action related to the 
reservoir plus interaction, 

S,l..«1=/*'{i(5-f -«--f ) (3.15) 

with 

00 

Hr^YI ^^n^nOin (3.16) 
n=l 

oo 
n,m=l 

In our case the lagrangean formalism simplifies the problem transforming a nonlinear 
problem into a linear one. The action (3.15) is quadratic in a, so it can be solved exactly. 
Observe that the Euler-Lagrange equations for (3.15) are 

oo 

din + i^nOI-n - Y g^nOtm ^ ^ (3.18) 
m=l 

oo 

a*^-iQ:nOin + ixY gnmOi*^ = ^ (3.19) 
m=l 

which must be solved subject to the boundary conditions 

an(0) = Pn (3.20) 

<it) = <■ (3.21) 
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Due to (2.17) we have gnn — and the modes are not coupled among themselves. This 
makes (3.18) and (3.19) easy to solve. That set of equations represents a set of harmonic 
oscillators forced by the presence of the soliton. The result can be written as, 

a„(r) = e-*""- + f; Wnm{r)(3m] (3.22) 

\ m=l / 



oo 



<(t) = e^^"- <e-^-* + ^ W^„^(T)e-'"-*< (3.23) 

\ m=l / 

where Wnm and Wnm are functionals of x{t) which obey the following equations 

CO ..^ 

Wr,^{T) = 5nm+Y. dt' Kr^At'Wn'mit') (3.24) 

^ ft 

Wnmir) ^Snm+T. / dt' K,,^>{t')Wn'm{t') (3.25) 
n—l 

where 

Knm{[x],r) = i5„^i(i')e^^""-'''"^*' (3.26) 

is the kernel of the integral equation (observe that: Wnmit) — M^mn(O)). 

Now we expand the action (3.15) around the classical solution (3.22) and (3.23) and 
obtain, after some integrations in (3.14) 

oo 

F[x,y] = ll{l-T^^[x,y]nn)-' (3.27) 



n=l 



where 



oo 

rnm = W:jy] + Wmn[x] + E WL[y]Wen[x] (3.28) 



with 



Notice that (3.27) and (3.28) are exact, no approximations have been made so far. 
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We see from (3.24) and (3.26) that Wnm can be expressed as a power series expansion of 
the Fourier transform of the sohton velocity, x, so for small soliton velocities we expect that 
only few terms in (3.24) will be sufficient for a good description of the soliton dynamics. 

Another way to see this is to notice that (3.24) and (3.25) are the scattering amplitudes 
from the mode k to the mode j. The terms that appear in the sum represent the virtual 
transitions between these two modes. With these two arguments in mind we will make use 
of the Born- approximation. In matrix notation, 



Therefore, in the approximation of small soliton velocity the terms in (3.28) are small 
and we can rewrite as a good approximation 



Observe that if the interaction is turned off (F 0) or the temperature is zero (T — 0) 
the functional (3.31) is one, and, as we would expect the soliton moves as a free particle. 
Substituting the Born approximation (3.30) in (3.31) and the latter in (3.12) we find 




(3.30) 




(3.31) 




(3.32) 



where 



jyt'[^{x\t')-fit')) + ixit')-yit')) [dr Vj{t'-r){x{r) + y{r))] (3.33) 



and 




(3.34) 



with 



00 



(3.35) 



n,m=l 



00 



(3.36) 



n,m=l 
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Now, if we define the new variables R and r as, 



R = ^ (3.37) 



r = X — y (3.38) 

the equations of motion for the action in (3.33) read, 

R{t) + 2 /* dt'-i{T - t')R{t') = (3.39) 
Jo 

f(r) - 2 /* dt'j{t' - r)r(t') = (3.40) 
Jo 

where 



or, using (3.36), 



7W = -r^ E ^?L^n(^^n-^^m)cos(l]„-l]Jt (3.42) 

^'^o n,m=l 



is the damping function. 

In terms of these newly defined variables, we can easily see that (3.39) and (3.40) have 
the same form of the equations previously obtained in the case of quantum brownian motion 
0, except for the fact that they now present memory effects (see ref. for details). 

Furthermore, in the limit where the time scale of interest is much greater than the 
correlation time of the meson variables we can write , 



7(t) = 7{T)Sit) (3.43) 

where 7(T) is a damping parameter which is temperature dependent and 6{t) is the Dirac 
delta function. The form (3.43) is known as the Markovian approximation, because in this 
case the memory is purely local and does not depend on the previous motion of the soliton. 

If we use (3.39) and (3.40) with (3.43) and expand the phase of (3.32) around this 
classical solution we get the well-known result for the quantum brownian motion P] where 
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the damping parameter 7 (temperature independent) is replaced by 7(T) and the diffusive 
part is replaced by (3.34). As a consequence, the diffusion parameter in momentum space 
will be given by, 

^ W = = -^'^('^) E almnni^n ' ^mf COs(fi„ - fi^)t (3.44) 

n,m=l 

As discussed before, in the Markovian limit, we can write D[t) in the Markovian form, 

D{t) = D{T)S{t) (3.45) 
where D{T) and 7(T) obey the classical fluctuation-dissipation theorem at low temperatures 



In what follows we shall define a function S{uj,uj') which will, in analogy to the spectral 
function J{uj) of the standard model 0, allows to replace all the summations over k by 
integrals over frequencies, 

oo 

Siu, u') = glJ{uj - n^)Sico' - ^]„). (3.46) 

n,m=l 

Notice, however, that unlike J{uj) in [Q, this new function S{uj,uj') is related to the 
scattering of the environmental excitations between states of frequencies u and u' (as seeing 
from the laboratory frame). Moreover, due to (3.24) it is easy to see that, 

S{u},uj') = S{uj',uj) (3.47) 

We call S{uj,Lj') the "scattering function". 

Notice that we can rewrite (3.42) and (3.44) as, 

7(t) = r du r duj'S{uj, uj'){uj - uj'){n{uj) - n{uj')) cos(cu - uj')t (3.48) 



2Mo Jo Jo 

and 
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oo roo 

r\2/ 



"'0 



D{t) = / duj duj'S{uj — uj '){uj - JY{n{uj) + n{J)) cos(u; - J)t. (3.49) 



Concluding, we have established that the Hamiltonian (2.27) leads to a brownian dy- 
namics, that is, the soliton moves as a particle in a viscous environment where its relaxation 
and diffusion are due to the scattering of mesons. 
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IV. CONCLUSIONS 



In this article we have showed how a new model for the study of dissipation in quantum 
mechanics can be obtained from the study the solitons in quantum field theory. We showed 
how the transport properties of solitons can be calculated in a consistent way from the 
microscopic point of view. 

Our results show that solitons move as a brownian particle at low energies due to the 
scattering of the mesons which are present in the system. These mesons are the residual 
excitations created by the presence of the soliton. We showed that the damping and the 
diffusion coefficients for the soliton motion are dependent of the temperature since the 
mesons must be thermally activated in order to scatter off the soliton. Therefore, at zero 
temperature the soliton moves freely, but its mobility decreases as the temperature increases 
(we have showed this explicitly for the case of the polaron motion, see ref. 0]). All these 
results can not be obtained from the standard model since the origin of the particle and the 
environment is distinct is that case. 
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